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APPLICATION  OF  PASTERNAK  MODEL  TO  SOME  SOIL-STRUCTURE  INTERACTION  PROBLEMS 


DETERMINATION  OF  FOUNDATION  PARAMETERS  AND  THE  CHOICE  OF  MODELS* 

PART  Is  INTRODUCTION 


1.  When  foundation  models  are  used  for  the  analysis  of  continuously 
supported  structures  that  occur  in  engineering  practice,  a  major  difficulty 
encountered  is  the  determination  of  the  foundation  parameter(s)  which  enter 
these  analyses. 

2.  The  four  main  methods  for  determining  these  parameters  found  in  the 
literature  are: 

Method  I  The  analytical  results  for  the  structure,  like  deflections 
or  bending  moments,  that  are  based  on  the  adopted  simple 
foundation  siodel,  are  compared  with  corresponding  test 
data  (References  1-4). 

Method  II  The  analytical  results  that  are  based  on  the  adopted  foun¬ 
dation  model  are  compared  with  the  corresponding  analyti¬ 
cal  results  for  which  the  foundation  model  is  an  elastic 
continuum  (References  5-9). 

Method  III  In  situ  plate  loading  tests  are  used  (References  10-17). 

Method  IV  The  foundation  response  expressions  and,  hence,  the  pa¬ 
rameters  are  derived  directly  from  the  equations  for  an 
elastic  continuum  by  introducing  a  priori  simplifying  as¬ 
sumptions  regarding  stresses  and  displacements  (Refer¬ 
ences  18-27)  or  by  a  formal  power  series  expansion  and 
retention  of  lower  order  terms  (References  20,  28-30). 

3.  The  fourth  method  has  gained  in  popularity  in  recent  years,  since  it 
suggests  that  the  parameters  of  a  simplified  foundation  model  can  be  calcu¬ 
lated  from  the  constants  E,  n  of  an  "equivalent"  elastic  continuum  which  are 
determined  from  laboratory  tests  on  small  soil  samples  taken  from  the  layer(s) 
of  the  subgrade  under  consideration. 


*  Research  supported  by  the  US  Army  Engineer  Waterways  Experiment  Station, 
Vicksburg,  Miss.  39180. 


4.  One  aim  of  this  paper  is  to  show  that  the  fourth  method  is  concep¬ 
tually  of  questionable  validity  and,  thus,  it  should  not  be  considered  for  the 
determination  of  foundation  model  parameters.  The  principal  aim  is  to  discuss 
the  criteria  for  choosing  suitable  foundation  models  and  the  methods  for 
determining  the  corresponding  foundation  parameters. 
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PART  II:  DERIVATION  OF  FOUNDATION  MODEL  PARAMETERS  AS 
FUNCTIONS  OF  THE  CONSTANTS  E,  v  OF  AN  ELASTIC  LAYER 


5.  According  to  Method  IV,  the  derivation  of  simplified  foundation  re¬ 
sponse  expressions  from  the  equations  of  a  continuum,  by  introducing  a  priori 
simplifying  assumptions,  directly  yields  the  coefficients  for  the  foundation 
response  expressions  in  terms  of  the  material  parameters  E,  V  and  the  layer 
thickness  H  .  These  foundation  response  equations  can  be  derived  by  any  of 
the  following  three  procedures: 

Procedure  1  Using  the  differential  equations  for  a  continuum  of 
finite  depth,  as  done  in  References  18-22,  or 

Procedure  2  Using  the  variational  approach,  as  done  in  Refer¬ 
ences  23-27,  or 

Procedure  3  Using  the  formal  approach  by  expanding  in  power  series 
and  retaining  terms  of  lower  order,  as  discussed  in 
References  20  and  28-30. 

6.  For  the  following  presentation,  Method  IV(l),  which  uses  the  differ¬ 
ential  equations  for  a  continuum  of  finite  depth,  will  be  used  because  the 

made  assumptions  have  a  simple  physical  meaning.  At  first  the  relevant  re- 

18  19 

suits  obtained  by  E.  Reissner  ’  are  summarized.  The  equations  which  corre¬ 
spond  to  the  Pasternak  foundation  and  to  the  Winkler  foundation  are  then 
derived.  This  is  followed  by  a  critique  of  the  determination  of  the  model 
parameters  from  the  derived  foundation  response  equations. 

The  Reissner  Foundation  Response 


7.  Consider  an  elastic  layer  of  thickness  H  which  rests  on  a  rigid 
base  and  is  subjected  to  a  distributed  load  p  ,  as  shown  in  Figure  1.  Assum¬ 
ing  the  simplifying  assumptions 


«„(**,*)  =  OyyCx.y.z)  =  a  (x,y,z)  =  0 

u(*»y»0)  =  v(x,y,0)  =  0 

u(x,y,H)  =  v(x,y,H)  =  w(x,y,H)  =  0 


(1) 


5 


-•  »-•  •  •• . 


Figure  1.  Elastic  layer  resting  on  1  rigid  base 

where  cr„  are  the  stresses  and  (u,v,w)  are  the  components  of  the  displace¬ 
ment  vector  in  the  direction  of  (x,y,z)  ,  Reissner  derived  the  foundation 
response  expression 

p  .  *d§  =  £  «  -  82  fy  C2) 

P12E  PH  3  K 

In  the  above  equation  W(x,y)  =  w(x,y,0)  is  the  vertical  deflection  of  the 

free  foundation  surface  and  p(x,y)  =  -  ct  (x,y,0)  is  the  distributed  pres- 

zz 

sure  normal  to  the  free  surface.  A  detailed  derivation  of  Equation  2  was 

22 

recently  presented  by  Horvath 

Derivation  of  the  Pasternak  Foundation  Response 
from  the  Equations  of  an  Elastic  Layer 


8.  Consider  the  elastic  layer  shown  in  Figure  1  and  assume  that  instead 
of  Equation  1,  the  simplifying  assumptions  are 


a  (x.y.z)  =  ct  (x,y,z)  =  a  (x,y,z) 
xx  yy  yx 


=  0 


(3) 


u(x,y,z)  =  0 


v(x,y,z)  s  0 


6 


xt,  we  substitute  these  assumptions  into  the  equilibrium  equations,  the 
ress-strain  equations,  and  the  strain-displacement  relations  for  a  linearly 
astic  solid  and  obtain 


9a 


9a  9  a 
9x  9y 


e  =  -  ^  a 
xx  E  zz 


-Kz 


e  =  |  a 
zz  E  zz 


e  =0 
xy 


2G  CTxz 


£yz  2G  °yz 


e  =0 
xx 


e  =0 

yy 


e  =0 
xy 


'xz  2  9z 


9.  The  first  two  equations  in  Equations  5  and  6  imply  that  v  has  to 
•  assumed  equal  to  zero  and  from  the  first  two  equations  in  Equation  U  it 
>llows  that,  as  in  the  Reissner  foundation,  the  shearing  stresses  are  inde- 
:ndent  of  z  .  Thus, 


a  =  a  (x,y) 
xz  xz 


a  =  a  (x,y) 
yz  yz 


i  an  unrealistic  result,  especially  for  thick  foundation  layers.  However, 
Lnce  foundation  models  are  usually  introduced  to  study  the  response  of  the 
>undation  surface  to  applied  loads,  and  not  the  stresses  and  displacements 
Lthin  the  foundation  material,  this  deficiency  may,  in  general,  be  of  no 
*rious  consequence. 

10.  Integrating  the  third  part  of  Equation  U  with  respect  to  z  ,  noting 
piation  7,  we  obtain 

/do  do  \ 

022<x’y’z)  =  -z  VTa?  +  “9^/  +  f(x’y)  (8) 


.  .N  > 


rroa  the  boundary  condition 


Lt  follows  that 


ind  thus 


azz(x,y,0)  =  -p  (x,y) 


f(x,y)  =  -p  (x,y) 


ozz(x,y,2)  =  -p  (x,y)  -  z 


i do  do  \ 

V  9x  3y  / 


11.  Eliminating  £ ^  from  the  third  equation  in  Equations  5  and  6,  we 

>btain 


3w  _  1 
dz  E  °zz 


Substituting  Equation  11  into  Equation  12,  and  then  integrating  the  resulting 
equation  with  respect  to  z  ,  we  obtain,  noting  Equation  7, 

2  (da  3o  \ 

E  w(x,y,z)  =  -z  p(x,y)  -  +  8<x»y)  03) 


From  the  boundary  condition 


w(x,y,H)  =  0 


it  follows  that 


„2  (da  do  \ 
g(x,y)  =  H  P(x,y)  ♦  f-  ) 


end  hence 


E  w(x,y,z)  =  (H  -  z)  p  (x,y)  + 


2  2 
H  -  z 


(da  da  \ 
__zx  +  _2y 

\  3x  ay  / 


According  to  Equations  5  and  6 


.  3w 
o  =  G  — 
zx  dz 


;  CT  =  G  p  and 
zy  3y 
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completed  vertical  walls,  and  then  to  analyze  the  walls  and  the  continuously 
supported  floor  plate. 

41.  In  the  past,  the  Winkler  foundation  model  was  used  to  represent  the 
response  of  the  base.  A  major  shortcoming  of  this  model  is  that  it  does  not 
represent  correctly  the  vertical  pressures  along  the  edges  and  (T)  .  Thus, 
the  resulting  bending  moments  in  the  floor  plate  may  deviate  strongly  from  the 
actual  ones. 

42.  To  eliminate  this  shortcoming,  either  the  Pasternak  foundation  or 
the  Kerr  foundation  may  be  used.  The  corresponding  models  of  the  structure  to 
be  analyzed  are  shown  in  Figures  8  and  9. 


43.  The  above  examples  demonstrate  that  the  choice  of  a  foundation 
model  depends  on  a  variety  of  factors  such  as  the  geometry  of  the  structure 
and  the  characteristics  of  the  subgrade.  In  this  connection  it  should  be 
noted  that  the  determination  of  foundation  parameters  is  meaningless  unless 
the  chosen  foundation  model  can  represent  the  anticipated  base  response  of 
interest. 
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Figure  4.  Model  for  railway  track  analysis 


Figure  6.  Plate  on  Kerr  foundation 


Figure  3.  Contact  pressure  distribution 


response  (for  example,  contact  pressure).  This  selection  should  be  based  on 
the  established  soil  profiles  of  the  location  under  consideration  and  the  type 
of  structure  to  be  built.  The  chosen  foundation  model  should  incorporate  the 
main  feature  of  the  anticipated  response. 

39.  For  example,  the  Winkler  foundation  represents  the  base  response 
rather  well  for  the  analysis  of  railway  tracks  subjected  to  the  static  wheel 
loads  of  locomotives  and  cars,  as  shown  in  Figure  4,  because  the  structure  is 
very  long.  It  is  not  suitable  for  representing  the  base  response  for  a  foot¬ 
ing  on  a  base  with  cohesion,  as  shown  in  Figure  3.  In  this  case  a  more  suit- 

31  32 

able  foundation  model  is  the  one  by  Pasternak  or  Kerr,  as  shown  in  Fig¬ 
ures  5  and  6. 

40.  In  designing  dock  structures  and  navigation  locks,  structures  of 
the  type  shown  in  Figure  7  are  encountered.  One  approach  for  analyzing  such 
structures  is  to  estimate  the  largest  net  pressures  to  be  expected  on  the 
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be  used  with  caution,  when  applied  to  actual  engineering  problems. 

32.  The  foundation  parameters  calculated  using  Method  III,  which  is  based 
on  in  situ  plate  loading  tests,  are  useful  only  if  the  resulting  analyses  will 
predict  the  corresponding  results  anticipated  in  the  field  (like  contact  pres¬ 
sures)  with  a  reasonable  accuracy. 

33.  Method  I,  according  to  a  collocation  of  analytical  results  based  on 
a  chosen  foundation  model  with  the  corresponding  test  data,  is  the  most  suit¬ 
able  one.  However,  it  requires  the  existence  of  the  structure,  or  a  large- 
scale  model.  This  method  is  used  for  the  determination  of  the  modulus  of 
railway  tracks.^ 

34.  It  should  be  noted  that  the  determination  of  foundation  model 
parameters  by  analytical  and  tested  collocation  of  surface  deformations  out¬ 
side  the  loaded  region  is  of  questionable  validity. 

Remarks  on  the  Choice  of  a  Proper  Foundation  Model 

35.  The  response  equation  is  usually  well  known  for  the  analysis  of  a 
continuously  supported  structure  in  geotechnical  engineering.  However,  diffi¬ 
culties  are  encountered  when  choosing  the  response  equations  for  the  support¬ 
ing  base. 

36.  The  reason  for  this  is  that  the  structure  is  made  of  a  solid  mate¬ 
rial  (metal  or  reinforced  concrete)  and,  while  in  service,  it  is  expected  to 
respond  elastically  when  subjected  to  static  and  dynamic  loads.  On  the  other 
hand,  the  subgrade  is  generally  a  granular  material,  often  with  a  variable 
degree  of  cohesion,  and,  hence,  may  easily  undergo  nonelastic  deformations  for 
any  type  of  load. 

37.  To  indicate  the  wide  range  of  possible  responses  of  a  subgrade, 

Figure  3  shows  two  contact  pressure  distributions  between  a  footing  and  its 
base  as  observed  in  tests.  Note  the  large  difference  in  the  distributions 
when  the  base  is  a  cohesionless  sand  or  a  base  with  cohesion.  The  contact 
pressure  also  varies  substantially  when  the  footing  rests  on  top  of  the  sand 
base,  or  is  acting  a  few  meters  under  the  base  surface.  Note,  also,  that  the 
elastic  continuum  predicts  a  very  different  pressure  distribution  from  the  one 
observed  for  a  sand  base. 

38.  Thus,  in  accordance  with  Step  2,  when  choosing  a  model  for  the  sub¬ 
grade  it  is  essential  to  realize  a  priori  the  nature  of  the  anticipated 
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Its  task  is  to  replace  the  defined  physical  problem  of  interest  by  an  ideal¬ 
ized  (simplified)  analytical  model  and  then  to  derive  the  formulation  for  it. 
Many  of  the  serious  deficiencies  of  engineering  analyses  originate  in  this 
step.  The  response  equations  for  continuously  supported  beam-  or  plate- type 
structures  are  well  known.  The  problematics  appear  when  choosing  the  re¬ 
sponse  equations  (model)  for  the  supporting  medium. 

27.  At  this  point  it  should  be  noted  that  an  exact  solution  of  a  formu¬ 
lation  yields  the  response  of  the  chosen  model,  but  not  necessarily  the  re¬ 
sponse  of  the  actual  structure.  Therefore,  a  desirable  objective  is  to  devise 
an  analytical  model  for  the  structure  and  foundation  whose  response  of  inter¬ 
est  is  very  close  to  the  actual  situation,  and,  at  the  same  time,  whose  formula¬ 
tion  is  easily  solvable. 

28.  Step  3  consists  of  the  solution  of  the  formulated  problem  and  this 
is  a  mathematical  exercise.  For  use  in  engineering  practice,  it  is  irrelevant 
whether  the  results  are  obtained  in  closed  form  or  numerically  using  a  computer 
method,  if  the  numerical  solution  exhibits  the  desired  degree  of  accuracy. 

29.  The  purpose  of  Step  4  is  to  determine  the  parameters  that  enter 
the  analysis  and  to  validate  the  derived  formulation,  two  tasks  that  are  very 
closely  related.  They  are  usually  performed  by  collocating  the  obtained  ana¬ 
lytical  results  with  corresponding  test  data  (Method  I).  This  validation 
should  determine  whether  or  not  the  derived  formulation  represents  the  problem 
under  consideration  with  the  desired  degree  of  accuracy,  including  a  wide 
range  of  parameters  which  may  be  encountered  in  engineering  practice;  the 
wider  the  range  of  parameters  the  better  the  model. 

Remarks  on  the  Determination  of  Foundation  Parameters 


30.  In  Part  III  it  was  shown  that  Method  IV  is  not  suitable  for  the  de¬ 
termination  of  foundation  parameters  used  in  engineering  practice. 

31.  Regarding  Method  II,  it  should  be  noted  that  the  determination  of 
foundation  parameters  from  a  collocation  of  a  solution  based  on  a  simple 
foundation  model  with  the  corresponding  solution  based  on  an  elastic  continuum 
is  often  only  of  academic  interest.  This  is  so  because  there  are  many  situa¬ 
tions  in  geotechnical  engineering  for  which  the  elastic  continuum  is  not  suit¬ 
able  for  representing  the  base  response  (for  example,  footings  on  sand). 
Therefore,  the  foundation  parameters  determined  by  using  Method  II,  should  also 
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PART  IV:  PRINCIPLES  FOR  CHOOSING  FOUNDATION  MODELS  AND  METHODS 
FOR  DETERMINING  THE  NECESSARY  PARAMETERS 


General  Considerations 

24.  To  begin  this  discussion  let  us  first  review  the  general  procedure 
for  obtaining  analytical  solutions  to  physical  problems.  This  procedure  con¬ 
sists  essentially  of  four  distinct  steps,  which  are  shown  in  the  following 
tabulation: 


Step  1 


DEFINITION 
OF  PHYSICAL 
PROBLEM 
to  be  analyzed 

For  example: 
Determination 
of  contact 
stresses  and 
bending  moments 
in  a  footing. 


Step  2 


FORMULATION 
OF  THE 
PROBLEM 
(Mechanics) 

1.  Choice  of  an 
idealized  model 

2.  Analytical 
formulation  of 
model 

Result:  Equations 


Step  3 


SOLUTION  OF 
FORMULATION 

(Mathematics) 

1.  Analytical 
solution  or 

2.  Numerical 
solution 


Step  4 


DETERMINATION 
OF  PARAMETERS 
AND  VALIDATION 
by  comparing  with 
results  of  rele¬ 
vant  tests 


(If  agreement  is 
poor,  return  to 
Step  2.) 


25.  Step  1,  although  superficially  obvious,  is  a  very  essential  step 
and  certainly  not  a  trivial  one.  This  is  true  because  a  very  general  formula¬ 
tion,  containing  a  large  range  of  information,  is  usually  difficult  to  derive 
and  even  more  difficult  and  costly  to  solve.  Therefore,  in  engineering,  as  in 
the  various  branches  of  the  physical  sciences,  it  is  often  necessary  to  limit 
the  range  of  answers  of  interest  to  be  obtained  from  an  analysis,  in  order  to 
simplify  the  resulting  formulation  and  solution  (namely.  Steps  2  and  3).  In 
this  approach,  if  more  than  one  answer  for  the  same  structure  is  needed  it  may 
require  aiore  than  one  formulation.  As  an  example,  note  that  the  analysis  of  a 
structure  and  its  foundation  should  be  based  on  one  formulation,  which  con¬ 
tains  the  equations  for  the  structure  and  the  equations  for  the  foundation. 
However,  in  order  to  simplify  the  necessary  analysis,  often  this  problem  is 
separated  into  two  parts:  (a)  a  structure  that  rests  on  a  simple  foundation 
amdel  and  (b)  a  continuum  foundation  that  is  subjected  to  the  structural 
forces  determined  from  part  (a) . 

26.  Step  2  is  usually  the  most  crucial  and  intricate  of  the  four  steps. 


Demonstration  Indicates  Value  Differences  in  Original 
Continuum  and  Simplifying  Assumptions 


22.  To  demonstrate  this  point  consider  an  analogous  situation  of  a  thin 
elastic  plate  that  is  governed  by  the  differential  equation  for  the  deflec¬ 
tions  w(x,y), 


DV  w  =  q 


where 


12  (1  -  v“) 

is  the  bending  stiffness  of  the  plate,  h  is  the  thickness  of  the  plate,  and 
(E,  v)  are  the  elastic  constants  of  the  plate  material.  Assume  that,  at  some 
instant,  the  plate  becomes  penetrated  by  a  dense  net  of  randomly  distributed 
fine  cracks  (possibly,  due  to  shrinkage  of  the  concrete).  It  is  easy  to  rea¬ 
lize  that,  in  this  new  state,  the  "effective"  bending  stiffness  D  of  the 
cracked  plate  is  lower  than  the  original  value;  the  reduction  depending  on  the 
density  and  type  of  the  crack  pattern.  Since  the  plate  thickness  h  did  not 
change  due  to  the  fine  cracks,  it  follows  from  Equation  31  that  the  "effec¬ 
tive"  E  and  v  differ  from  the  original  values.  The  fine  crack  pattern  in 
the  plate  problem  is  analogous  to  the  simplifying  assumptions  made  when  de¬ 
riving  the  foundation  response  expressions  from  the  equations  of  an  elastic 
continuum.  Thus,  the  E  and  v  values  of  the  "continuum"  that  is  subjected 
to  the  simplifying  assumptions  are  not  the  same  as  those  of  the  original 
continuum . 

23.  The  above  discussion  suggests  that  the  procedure  for  determining 
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the  foundation  modulus  as  proposed  by  Horvath,  as  well  as  by  some  other 
researchers  who  used  Method  IV,  is  generally  not  suitable  for  determining  the 
foundation  parameters  for  use  in  engineering  practice. 


w,~t 


19.  Next,  consider  the  nonhomogeneous  foundation  discussed  by  Horvath/ 
naaely  when  E  =  a  +  pz  ;  a  end  p  being  constants.  Also,  for  this  case, 
differential  Equation  24  can  be  integrated  directly  to  yield 


,(2)  *  r/s-fu iz  * 


Setting  a  +  pz  =  £  and  nondiaensionalizing  the  denominator,  it  follows  that 


w<*)  *  p/t  d  (£) 


+  ®2  =  0A  in  +  +  B2 


The  integration  constants  and  are  determined  from  the  boundary  con¬ 

ditions  in  Equation  25.  The  resulting  displacement  expression  is 


/  x  _  P  ,  a  +  pH 

”(z)  s* Jn  5-rfe 


0  <  z  <  H 


At  z=0,  with  P/A  =  p  and  w  =  W  ,  Equation  28  can  be  written  as 


p  -  2n”(7«  +~0H T/oJ  9 

This  expression  is  the  saaie  as  Equation  27,  except  for  the  coefficient  of  W  . 
The  coefficient  of  W  is  identical  to  the  k  expression  given  in  Refer- 
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ence  21,  Equation  9a,  for  E  =  a  +  pz  ,  which  was  derived  froa  the  "simplified 
continuua  procedure." 

20.  The  above  derivations,  Equations  24-29,  demonstrate  that  E/H  and 
P/2n  ((a  +  pH) /a]  are  the  k-values  for  a  simple  coluan  of  height  H  .  Thus, 
each  is  the  foundation  aodulus  for  a  base  that  consists  of  closely  spaced 
independent  vertical  columns  (springs).  However,  this  does  not  iaply  that 
each  is  the  foundation  aodulus  k  for  a  truly  continuous  base,  as  encountered 
in  engineering  practice. 

21.  When  reviewing  attesq>ts  to  determine  foundation  parameters  from  the 
constants  E,  v  of  an  elastic  continuua,  as  indicsted  in  Equation  20  or  23, 
it  has  to  be  noted  that  the  introduction  of  simplified  assumptions,  of  the 
type' stated  in  Equation  1,  3,  or  21,  changes  the  nature  of  the  continuum,  and 
affects  the  constants  E  and  v  . 


Figure  2.  Foundation  layer  consisting  of  closely 
spaced  vertical  columns 


d_ 

dz 


(0  <  z  <  H) 


and  the  boundary  conditions 


w' (0)  =  -  (P/EA)z=q 
w(H)  =  0 


(2- 


(2! 


In  above  formulation  w(z)  is  the  vertical  displacement  at  z  ,  A  is  the 
cross-sectional  area  of  the  column,  and  E  is  Young's  modulus  of  the  column 
material. 


18.  For  the  case  when  EA  =  constant  ,  the  above  fonaulation  yields 


w(z) 


P(H  -  z) 
EA 


(2< 


At  z  =  0  ,  with  P/A  =  p  and  w  =  W  ,  Equation  26  can  be  written  as 


(2; 


which  is  identical  with  Equation  22. 


PART  III:  COMMENTS  TO  THE  DETERMINATION  OF  THE  MODEL  PARAMETERS 
FROM  THE  DERIVED  FOUNDATION  RESPONSE  EXPRESSIONS 


Horvath  Approach 

15.  The  parameter  expressions  of  the  type  given  in  Equations  20  and  23 

led  a  number  of  researchers  to  the  conclusion  that  the  foundation  model  param*- 

eters  may  be  calculated  from  the  constants  E,  v  of  an  equivalent  elastic 

continuum,  which  are  determined  from  laboratory  tests  on  samples  taken  from  the 

actual  base.  The  most  recent  example  of  this  approach  is  described  in  the 
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paper  by  Horvath.  He  derived  the  foundation  response  expression 

P(x,y)  =  |  W(x,y)  (22*) 

as  indicated  in  paragraph  12,  noted  that  the  above  equation  is  identical  to 
the  Winkler  foundation  response  when  the  Winkler  parameter  is 

k  =  |  (23  * ) 

and  then  suggested  “...that  k  can  be  interpreted  and  evaluated  using  an 
approach  that  has  as  its  basis  the  theory  of  elasticity..."  and  thus,  that 
“...it  is  possible  to  evaluate  k  from  known  soil  parameters..."  E,  v  and 
the  layer  thickness  H  ,  using  Equation  23. 

Procedures  Proposed  by  Horvath  and  Other 
Researchers  Conceptually  Invalid 

16.  It  is  the  opinion  of  this  writer  that  the  above  procedure  proposed 
21 

by  Horvath  is  conceptually  of  questionable  validity.  This  applies  also  to 
similar  suggestions  made  by  some  of  the  other  researchers  who  used  Method  IV. 
(This  objection  does  not  apply  to  the  determination  of  foundation  parameters 
in  terms  of  E,  v  used  in  Method  II,  which  is  based  on  a  different  concept.) 

17.  To  illustrate  this  point,  note  that  the  simplifying  assumptions  in 
Equation  21  suggest  that  the  foundation  layer  consists  of  closely  spaced  ver¬ 
tical  columns  of  length  H  ,  as  shown  in  Figure  2,  that  are  not  interconnected. 
The  standard  analytical  formulation  of  these  columns  consists  of  the  differ¬ 
ential  equation 


with  the  stipulation  that  v  =  0  ,  as  in  the  previous  section. 

14.  Proceeding  as  before,  we  obtain  for  the  free  surface  of  the  founda 
tion,  instead  of  Equation  19,  the  response  expression 

(22 

This  expression  is  identical  to  the  response  equation  of  the  Winkler  founda¬ 
tion,  if  we  set 


noting  above  relations,  Equation  16  assumes  the  form 


E  w(x,y,z)  =  (H  -  z)  p  (x,y)  +  §  (H2  -  z2)  ^w  (18) 

_2  2  2  2  2 

where  v  =  8  /8x  +  8  /8y  is  the  Laplace  operator. 

12.  At  the  free  surface  of  the  foundation  z  =  0  and  Equation  18  be¬ 
comes  ,  setting  w(x,y,0)  =  W(x,y) 


(19) 


p(*»y)  =  |  w  “  ^ 


This  expression  is 
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tion,  if  the  two 


identical  to  the  response  equation  of  the  Pasternak  founda-  | 

foundation  parameters  are  set  j 

I 


k 


P 


E 

=  H  5 


(20)  | 

I 
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Derivation  of  the  Winkler  Foundation  Response 
from  the  Equations  of  an  Elastic  Layer 

13.  Consider  again  the  elastic  layer  shown  in  Figure  1,  but  assume, 
instead  of  Equation  3,  the  simplifying  assumptions,  that  throughout  the 
foundation 

o  =  a  =  a  -a  =  a  =  0  I 

xx  yy  xy  xz  yz  f 

>  (21) 

u  =  v  =  0  1 

With  these  assumptions  the  equilibrium,  stress-strain,  and  strain-displacement 
relations  for  the  elastic  layer  reduce  to  the  following  necessary  equations: 


8a 


zz 

8z 


=  0 
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E  ~ 
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